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Abstract 

We present an accurate quantum mechanical study of molecule-molecule collisions in the presence 
of a magnetic field. The work focusses on the analysis of elastic scattering and spin relaxation 
in collisions of 02('^5]~) molecules at cold (~0.1 K) and ultracold ('^lO^^ K) temperatures. Our 
calculations show that magnetic spin relaxation in molecule-molecule collisions is extremely efficient 
except at magnetic fields below 1 mT. The rate constant for spin relaxation at T = 0.1 K and a 
magnetic field of 0.1 T is found to be as large as 6.1 x 10~^^ cm^/s. The magnetic field dependence of 
elastic and inelastic scattering cross sections at ultracold temperatures is dominated by a manifold 
of Fcshbach resonances with the density of ~100 resonances per Tesla for collisions of molecules 
in the absolute ground state. This suggests that the scattering length of ultracold molecules in 
the absolute ground state can be effectively tuned in a very wide range of magnetic fields. Our 
calculations demonstrate that the number and properties of the magnetic Feshbach resonances are 
dramatically different for molecules in the absolute ground and excited spin states. The density 
of Fcshbach resonances for molecule-molecule scattering in the low-field-seeking Zeeman state is 
reduced by a factor of 10. 

PACS numbers: 33.20.-t, 33.80.Ps 
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I. INTRODUCTION 



The experimental work with cold and ultracold molecules is predicted to lead to many 
fundamental applications in quantum computation [1], condensed- matter physics [2], preci- 
sion spectroscopy [3], and physical chemistry [4]. Recent theoretical work has shown that 
ultracold ensembles of molecules trapped in optical lattices can be used to develop novel 
schemes for quantum information processing [6], design experiments for quantum simula- 
tions of condensed matter physics [2, 7], engineer novel phases with topological order [8], 
and explore many-body dynamics of strongly interacting systems [9]. Polar molecules in 
external field traps may form chains, which can be used to study rheological phenomena 
with non-classical behavior [2]. The creation of a Bose- Einstein condensate of molecules 
may enable the study of Bose-enhanced chemistry [10] and the effects of symmetry breaking 
on chemical interactions at ultracold temperatures [11]. The realization of these proposals as 
well as the creation of dense ensembles of ultracold molecules depends critically on the pos- 
sibility of controlling binary molecule - molecule interactions in molecular gases by external 
electromagnetic fields. The discovery of magnetic Feshbach resonances in atomic coUisions 
[12, 13] opened the door to many groundbreaking experiments with ultracold gases such as 
the realization of the BEC-BCS crossover [14], the observation of quantum phase transitions 
[7] and the creation of ultracold molecules using time- varying magnetic fields [15]. Extension 
of this work to molecular collisions may similarly lead to the development of new research 
directions such as cold controlled chemistry [4], quantum coherent control of bimolecular 
reactions [16], and quantum condensed- matter physics with molecular condensates [5]. 

Fueled by the promise of new discoveries, the experimental research of cold and ultracold 
molecular collisions is expanding rapidly. Some recent landmark experiments include the 
observation of threshold collision laws in Xe-OH scattering using slow molecular beams [17], 
measurements of cross sections for D2-OH collisions in a magnetic trap [18], and the detec- 
tion of magnetic Feshbach resonances in collisions of Cs2 molecules [19]. Further progress 
in the experimental study of complex molecule-molecule collisions and chemical reactions 
at cold and ultracold temperatures requires rigorous theoretical calculations elucidating 
the mechanisms of energy transfer in ultracold molecular collisions. Theoretical studies 
of molecule-molecule collisions at low temperatures are also necessary to understand the 
prospects for evaporative cooling of molecules in a magnetic trap. The evaporative cooling 
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relies on elastic molecule - molecule collisions as the trap depth is gradually reduced [13]. In 
order to sustain efficient evaporative cooling down to quantum degeneracy, the ratio of the 
probabilities for elastic scattering and spin relaxation (7) in molecule - molecule collisions 
must exceed 10^ [13]. The magnitudes of the rate constants for elastic and inelastic collisions 
and the dependence of the scattering observables in molecule - molecule encounters on the 
magnetic field are, however, completely unknown. 

Previous theoretical work by Krems and Dalgarno [20] and Volpi and Bohn [24] identified 
the main mechanisms of spin relaxation in atom-molecule collisions. It was found that spin 
relaxation in collisions of molecules with He atoms occurs through coupling to rotationally 
excited states mediated by the spin-spin interaction [20]. While the mechanisms of spin 
relaxation in molecule-molecule collisions should be similar, an accurate computational study 
of molecule-molecule scattering in a magnetic field is urgently needed to elucidate the rates 
of energy relaxation in dense ensembles of trapped molecules and understand the prospects 
for evaporative cooling of molecules. The quantum dynamics of ultracold molecule-molecule 
collisions in the absence of external fields has been studied by several authors (see, e.g., 
Refs. [21-23]). Avdeenkov and Bohn calculated the rates of spin relaxation in collisions of 
02('^S~) molecules at zero magnetic field and observed Fcshbach resonances due to molecular 
rotation [25, 26]. Krems and Dalgarno [20] presented a formalism for quantum scattering 
calculations of cross sections for molecule-molecule coUisions in a magnetic field. However, 
they did not consider the symmetrization procedure required to properly describe collisions 
between identical molecules. Here, we extend the approach of Krems and Dalgarno to study 
collisions of identical molecules. The use of exchange symmetry reduces the number of 
scattering channels by a factor of two, which allowed us to obtain converged cross sections 
for elastic scattering and spin relaxation of 02(^S~) molecules in an external magnetic field. 

We consider ultracold collisions of O2 molecules in their ground electronic state of ^E~ 
symmetry. O2-O2 collisions play an important role in atmospheric chemistry [29, 30]. The 
O2-O2 dimer is an interesting molecular complex [27-35] (for a recent review, see Ref. [34]). 
The binding in the complex is affected by the Heisenberg exchange interaction and the 
formation of incipient chemical bond [31, 32]. Molecular beam scattering experiments [31, 32] 
and high-level ah initio calculations [33] have previously been used to elucidate the nature of 
the chemical bond in the O2-O2 dimer. Here, we calculate the magnetic field dependence of 
the probabihties for spin relaxation in O2-O2 collisions at both ultracold (down to 10~^ K) 
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and cold (0.5 K) temperatures relevant for magnetic trapping experiments. Our calculations 
elucidate the possibility of evaporative cooling of paramagnetic molecules in a magnetic 
trap and indicate that spin depolarization in molecular collisions can be controlled by an 
external magnetic field. The results of our work will be useful for the interpretation of 
collision experiments involving slow beams of O2 molecules produced by Zeeman deceleration 
[37] and cryogenic coohng [38]. 

II. THEORY 

A. Hamiltonian and symmetrized basis functions 

The Hamiltonian for two molecules in an external magnetic field can be written as 
[20, 25] 

where the vectors Va and re describe the orientation of the molecules A and B in the space- 
fixed frame, R is the vector joining the molecules, i is the orbital angular momentum of the 
collision complex, fi is the reduced mass, and V^i is the electrostatic interaction potential 
(see Sec. JIB). The asymptotic Hamiltonian describes the individual molecules A and 
B in the presence of a magnetic field [20]. 

-^as = Ha + Hb, (2) 
The Hamiltonians of an individual ^E~ molecule can be written as [39] 

= B,Nl + 2/xbS • + 7^. • + 3 Ass f ^ j ^s^l^-) [S. ® S,]f , (3) 

where the index v — A,B labels the molecule, is the rotational angular momentum, 
is the rotational constant, /ib is the Bohr magneton, B is the magnetic field vector and Si, is 
the electron spin. The spin-rotation and spin-spin interactions in Eq. (3) are parametrized 
by the phenomenological constants 7 and Ass, which do not depend on i/ since we consider 
collisions of identical molecules. The spin-spin interaction given by the last term in Eq. (3) 
can be represented as a sum of tensor products of two identical spin operators and spherical 
harmonics, which depend on the orientation fi, of the diatomic molecule in the space-fixed 
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(SF) coordinate frame [39] . We use the rigid rotor approximation — f„) and neglect the 
weak magnetic dipole-dipole interaction and the hyperfine interaction of ^^O [25]. 

We solve the scattering problem for two identical bosonic molecules by expanding the total 
wave function in a space- fixed uncoupled basis introduced in Ref. [20] . The uncoupled basis 
functions are products of the eigenfunctions of the operators and S*^, their projections 
along the space- fixed ^-axis N^,^ and S^^^ and the spherical harmonics |£m^) 

iTATB^m^) = |rA)|rB)|^m^) (4) 

where \ti,) — \Ni,Mj^^)\S,jMs^) ■ The matrix elements of the Hamiltonians (1), (2) in the 
basis (4) were derived by Krems and Dalgarno [20] for the case of distinguishable molecules. 
The basis functions for identical molecules must be the eigenfunctions of the permutation 
operator P : va Vb'iVb rA-,R^ — -R, which commutes with the total Hamiltonian 
(1). The symmetrized basis functions can be constructed by applying the operator 1 -|-P to 
Eq. (4) and normalizing the result [20, 22, 40, 41] 

where e = (— )^ and the symmetry of the basis functions with respect to the interchange 
of identical molecules is given by index rj which is equal to 1 for composite bosons and 
— 1 for composite fermions [42]. In what follows we assume that the first basis function in 
each product on the right-rand side of Eq. (5) depends on the coordinates of molecule A, 
and the second basis function depends on the coordinates of molecule B [43]. Equaton (5) 
defines the well-ordered (combined) molecular states with Ta > tb- The normalization factor 
[2(1 -|- ^tatb)]^'^^ takes into account the fact that the basis function (5) with ta = tb need 
not be symmetrized when rje — +1 (i.e. when the colliding molecules are in the same state). 
The symmetrized function (5) with ta = tb vanishes identically when rje — —1. 

The basis functions given by Eqs. (4) and (5) are also the eigenfunctions of the inversion 
operator, with eigenvalues given by (^—^^a+Nb+i |^4]^j gin^g fop homonuclear molecules both 
Na and A^b ai'e cither even or odd, the quantity e = (— )^ is conserved as well. Thus, it is 
convenient to rewrite Eq. (5) in a factorized form 

<PZrBern,=^'rlr,(rA,rB)\irne) (6) 

where 

^Ir^ {ta, re) = [2(1 + &^^^)] 1/2 O^aTb) + r?e|TBTA)] . (7) 

6 



We note that no such symmetry exists for heteronuclear molecules, where the basis functions 
with A£ = ±1 are coupled by the interaction potential, and e is not conserved. 

The uncoupled symmetrized basis functions given by Eq. (5) are not the eigenfunctions of 

the asymptotic Hamiltonian (2) because the spin-rotation and spin-spin interactions couple 
the states with different N, M]sr, and Ms- Since the scattering ^-matrix must be defined in 
terms of the eigenfunctions of Hgs, it is necessary to transform the wave function to a new 
basis [20] 

^l7B(^A,rB) = Yl C:?:^,7A7B^I^(^A,rB), (8) 

T-A>TB 

where the coefficients C'tItb.tatb ^^"1 the matrix, which diagonalizes the asymptotic Hamil- 
tonian in the symmetrized basis (5). More explicitly, C^HagC = E, where E is the diagonal 
matrix of asymptotic energies and e-y^-y^ = e-y^ -|- e-yg . The energies of individual molecules 
e^y can be obtained by diagonalization of individual molecule Hamiltonians H^, (2) in the 
subspace of functions \ti,). Exploiting the properties of the interchange operator P"^ = P 
and P^ — 1 and using Eq. (5), we obtain the following expression for the matrix elements 



of 



{Kr^ime\Hs.\(l>T^ri,e'm') ^ [(^ ^ ^, )]l/2 ['^re.r^ (ta| /^Alri) +7^6(^,3, (taI/ZaIt^)] ■ 

(9) 

This expression provides a convenient method of constructing the matrix elements of the 
asymptotic Hamiltonian (2) in the symmetrized basis. Alternatively, one can obtain the 
transformation coefficients (8) directly from the eigenvectors of Ha defined as 

|7a) = J^a^TAM. (10) 

The transformation can be derived by multiplying Eq. (10) by a similar expression for |7b) 
and rearranging the terms. The result is 

^rArB,7A7B ~ TTTTlA VTTlA V|T72 [^'^A7A^-rB7B + ^^^rA7B^TB7A] (H) 



[(l + 5rA,rB)(l + <57A. 



IB J 



We have verified by numerical tests that Eqs. (8) and (11) give identical results up to an 
unimportant overall phase. In practice, however, it is more convenient to use Eq. (11) 
because it requires the evaluation and the diagonalization of smaller matrices. 
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B. Interaction Potential 



The spin-dependent interaction potential between two molecules can be represented 
in the form 

1 s 

V4R,fA,rB) = Yl J2 ^s{R,rA,rB)\SMs){SMs\, (12) 

5=-l Ms=-S 

where S = Sa + Sb is the total spin and Ms — Msj^ + Ms^ ■ The representation (12) is often 
used to describe the interactions of ultracold alkali metal atoms [12]. Here, we choose an 
alternative approach proposed by van der Avoird and coworkers [27, 28], in which the full 
interaction potential is separated into the spin-independent and spin-dependent parts 

Vei{R, rA, fe) = Ki{R, rA, tb) + V;d(-R, ^a, tb). (13) 

The spin-dependent part, also known as the Heisenberg exchange interaction [27, 28], can 
be parametrized for two molecules as 

Vsd{R: rA, re) = -2J{R, ta, fB)^A • ^b, (14) 

where J{R,rA,rB) is a scalar function similar to the interaction potential but decreasing 
exponentially with R [28]. We expand the interaction potential in the angular basis [20, 27, 
44] 

^si(i«,rA,rB) = (47r)3/2 J] Vx^xMAx^Xsx{R,fA,f-B), (15) 

Xa,Xb,X 

with the space-fixed angular basis functions defined as 

^AAABA(^,rA,rB) = I \Yxj,m>,JfA)Yx^m^JfB)Yxm^{R), (16) 

■mxj^,mx^,mx \'^Xa ^Ab ^A J 

where (:::) denotes a 3-j symbol. An efficient procedure for the evaluation of the radial 
expansion coefficients in Eq. (15) is described in Appendix A. 

The matrix elements of the spin-indepenent interaction potential between the sym- 
metrized basis functions (5) have the form 

(CrB^mjHil'/'rkrm^) = [(1 + 5^^^) (1 + ^, )] 1/2 (1^) 

X [{TArB£m\Vsi\T'jyT^e'm'^) + r]e{TATBimi\Vsi\T^T'^£'m'^)] 
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The first (direct) term on the right-hand side is a matrix element in the unsymmetrized 
basis. Krems and Dalgarno [20] evaluated this matrix element using the expansion of Vgi 
in terms of the angular basis functions (15). Using the conservation of the total angular 
momentum projection M = Mj^^ + ^s'a '^^Nb +-^Sb i^i collisions in a magnetic field, 
we obtain 

V^A.AeA(i?)[(2WA + l)(2iV; + l) 

Aa,Ab,A 

X (27Vb + l)(2iV^ + 1)(2£ + 1)(2£' + 1)(2Aa + 1)(2Ab + 1)(2A + i)]V2(_)M^^+M^3+m, 

Aa Ab a \ / TVa Aa W 

Mn^ - ^atb - M'j,^ me-m'J \ -Mj,^ M^^ - M'j,^ M'^J \ Q 

iVB Ab iV^ W iVe Ab iV^ W £ A f W £ A f 

-M^3 M^Tb - Mjv3 Mjv3 y y y y -m^ - y y 

(18) 

The second (exchange) matrix element appears as a result of the symmetrization procedure 
(5). It can be obtained from Eq. (18) by interchanging the indexes A^a ^ ^'na ^ -^a^b' 
and M's^ ^ Mg^. 

The matrix elements of the Heisenberg exchange interaction are obtained using the same 
procedure: 

{^ATBtrn^KMlr^E'ru') = [(1 + ^^^^ ) (1 + ^, )] 1/2 (l^) 

X [{Sp^MsA{SbMsb\Sk ■ SB\SKM's^)\SBM'sJ{Tp^TBlm,\ - 2JAB|TX/m^) 
+ r]e{SAMsA{SBMsB\SA ■ -5B|-SBM^J|5AM^^)(TATB^m,| - 2 JABlr^ri/m^)] . 

The matrix elements of the prefactor —2J{R, 9a, Ob, <p) are given by Eq. (18) with VaaAba(-R) 
replaced by —2Jxj^\j^\{R)- Equation (19) shows that unlike the spin-independent matrix 
elements (17), both the direct and exchange contributions to the Heisenberg exchange matrix 
element are not diagonal in spin quantum numbers Ms^ and Ms^^. For completeness, we 
give here the expressions for the matrix elements of the operator Sa- Sb- The direct matrix 
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element has the form [20] 

{SaMs^\{SbMsJSa ■ Sb\SaM'sJ\SbM'sJ = 5ms^m',Jms^m',^Ms^Ms^ 
+ ISms,,m's,±iSms,,m's,^i [Sa{Sa + 1) - M's^iM's, ± 1)] [Sb{Sb + 1) - M^^ (M^^ t 1)] 

(20) 

and the exhange matrix element can be obtained from Eq. (20) by interchanging the indexes 
M's^ and M's^. 

C. Close-coupling equations, scattering amplitudes, and cross sections 

The symmetrized wave function of the collision complex can be expanded as 

= E E Fl,,ernSmrne{R)^,:,,{rA, Tb). (21) 

7A>7B t,mi 

Substituting this expansion into the Schrodinger equation, Hijj^ = Ei/j'^, where E is the total 
energy, we obtain the close-coupling equations for the radial expansion coefficients 

r eje + i) 

The asymptotic form of the solutions F^j^^^emtY 7' i'm' (-^) R ^ ^ defines the symmetrized 
scattering matrix 

^lAlYiimtaWm'.iR ^ 00) ^ Kl'A^l^l'Ju'^m,m', e^vl-^k^Al^R - 7r£/2)] 

- irr^) '''^^A.e.n...;.^..^ exp[i(^,, - .e/2)] (23) 
The scattering amplitude for indistinguishable molecules can be written as [23, 52] 

where the scattering amplitudes on the right-hand side can be written in terms of unsym- 
metrized T-matrix elements 

57A7B^7A7B(-^i' ^ 27r ^ ^ Yg^^(A)Yg/^^(^)r^^^3fo„^.^^^^£/„^. (25) 



^lAjBtme (-^) - E E [C^^C] 7ATB^m^ ;7i7^^^'m^-^7i7^5^'mJ ' ^^^^ 
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Substituting this expression into Eq. (24) and rearranging the terms, we find 



where 



^7A7B^mf;7^7^j^'m^ ~ ^(^l ^ ^(^l ^, ^]l/2 [^7A7B^m^;7A7B^'"^f + ^^^7A7B^m^ ;7B7A^'f"J " C^'^) 



1 

These T-matrix elements are exactly the same as those obtained from the solution of the 
close-coupling equations in the symmetrized basis set (as described above). They are related 
to the yS'-matrix elements via [48] 

^7A7B^m^;7A7B^'K ^ '^7A7^v'^7B7B'^^^"^'"^"»f ~ '^7A7B^m^;7A7B^'m^ " ^^^^ 

The integral cross section for indistinguishable molecules in fully polarized nuclear spin 
states can be obtained from the scattering amplitude as described in Appendix B. The final 
result is 



tt(1 + 6. 



7A7b; y^y^ l^r; |2 (2q\ 



'7A7b-^7a7b — p I 7A7B^m^;7A7B^H 

TA7B e^jne e'm'^ 

where /c^^^g = 2//(£' — e^^ — e^g) is the wave vector for the initial collision channel. This 
equation shows that the elastic cross section of indistinguishable bosons (or fermions) is 
twice as large as that of distinguishable particles. Our Eq. (29) agrees with the results 
presented by Avdeenkov and Bohn [25] and Burke [47]. 



D. Computational details 

To parametrize the asymptotic Hamiltonian (2), we used the accurate spectroscopic con- 
stants of ^'^02 (in units of cm"^): Be = 1.353, 7 = -0.00396, Ass = 1-985 [53]. In order to 
obtain converged cross sections at collision energies below 0.5 K, we used a basis set compris- 
ing 3 rotational states {Na, A^b = — 4) and 4 partial waves (£ = — 6). These parameters 
resulted in 2526 coupled differential equations for M = 0. At the lowest collision energy of 
10~^ K, four partial waves were sufficient to achieve convergence. The coupled equations 
(22) were solved on a grid of R from 2 ao to 200 uq with a step size of 0.04 qq yielding the 
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cross sections converged to within 10%. The spin- independent and spin-dependent parts of 
the interaction potential for the oxygen dimer were constructed to reproduce the integral 
cross sections measured in molecular beam scattering experiments [31, 32]. The lowest- 
order isotropic part Vooo(-R) was determined with high accuracy from an analysis of the 
glory pattern in high- velocity collisions of rotationally "hot" O2 molecules. The anisotropic 
coefficients V202(-R), V22o{R), and V222(-R) were inferred from the scattering experiments with 
rotationally cold, supersonic beams of aligned O2 [31, 32]. Recent ab initio studies [35, 36] 
have shown that the O2-O2 interaction may be more anisotropic than suggested in Ref. 
[31]. The inelastic transition probabilities in molecule-molecule collisions increase with the 
anisotropy of the interaction potential [20] , so the cross sections for spin relaxation presented 
in this work should be viewed as lower bounds to the actual magnitudes. 



III. RESULTS AND DISCUSSION 



We consider coUisions of two identical oxygen molecules in the ground electronic state 
^E~. We choose the ^^O isotope because the ^^02 molecule is characterized by the Ja = 1 
ground state required for magnetic trapping (where Ja = iVA + Sa is the total angular 
momentum of O2 exclusive of nuclear spin). The Zeeman levels of ^''©2 obtained by diago- 
nalizing the asymptotic Hamiltonian are shown in Fig. 1(a). Magnetic fields split the J\ = 1 
ground state of ^^62 into three Zeeman states corresponding to M^^ — 0, —1, and -|-1. We 
first consider coUisions of two oxygen molecules in the spin-aligned state \ Msj^ — 1, Msj^ — 1) , 
relevant for magnetic trapping and cryogenic cooling experiments. Molecules prepared in 
this state can scatter elastically or undergo spin-changing inelastic transitions. The lower 
panel of Fig. 2 shows that spin-changing collisions are suppressed at low magnetic fields (< 1 
mT) leading to a large value of the elastic to inelastic ratio 7 ~ 10"^ which is favorable for 
evaporative cooling. Since M is conserved, the spin relaxation transitions must be accom- 
panied by the transition i — ^ i' — 2, which leads to a centrifugal barrier in the outgoing 
reaction channel suppressing the s-wave inelastic scattering. Because the i — 2 barrier has 
a fixed height of 13 mK [25], this suppression occurs only at low magnetic fields where the 
energy difference between the initial and final Zeeman states does not exceed the barrier 
height [24, 25]. At larger magnetic fields, the spin-changing transitions are much more effi- 
cient, as shown in Fig. 2(6). This suggests that evaporative cooling of molecules in a deep 
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magnetic trap will be challenging. However, our results indicate that evaporative cooling 
may be possible in shallow magnetic traps (< 1 mT deep). The maximum temperature of 
the molecules that can be held in such a trap is determined by the parameter 

where T is the temperature, k-Q is the Boltzmann's constant, and // = 2^q for molecules. 
Efficient trapping requires 77 > 5 [59], which corresponds to a temperature of 0.25 mK for 
molecules aX B — 1 mT. Such temperatures can be reached using optical Stark deceleration 
[58] or Zeeman slowing [37]. 

Spin relaxation in collisions of S-state molecules in fully spin-stretched states can only 
occur through coupling to the rotationally excited states assisted by the spin-spin interaction 
[20]. This is a two-step mechanism illustrated schematically in Fig. 1(6). The spin-spin 
interaction (denoted by Vss in Fig. 1) mixes the ground rotational state = 0, Mat = 
0,M5 = 1) with the excited rotational states \N — 2,Mn,Ms — —1). Spin relaxation is 
then induced by the anisotropy of the interaction potential leading to the \N = 2, M^, Ms = 
— 1) \N' = 0, M'j^, M'g = —1) transition. The anisotropy of molecule - molecule interaction 
potentials is usually very strong; our results yield the rate constant for spin relaxation 
6.1 X 10-1^ cmVs at T = 0.1 K and S = 0.1 T. 

Figure 3 shows the magnetic field dependence of O2 - O2 scattering at a coUision energy of 
10~^ K. The spin relaxation cross section increases from zero to a large value over a narrow 
region of magnetic fields 5 = — 1 mT. For magnetic fields larger than 10 mT, the spin 
relaxation cross section summed over all final Zeeman states is a factor of 10 larger than the 
elastic scattering cross section. The cross section displays several broad resonances. Figure 
1(6) illustrates that at certain magnetic fields, the energy in the incoming collision channel 
becomes degenerate with that of a quasibound state supported by the uppermost curve. 
This results in the formation of a long-lived complex in which one of the molecules is in 
the N = 2 rotationally excited state [25] . The lifetime of the complex is determined by the 
strength of the couplings induced by the spin-rotation interaction and the intermolecular 
potential as well as the presence of inelastic loss processes. In particular, the peaks shown 
in Fig. 3 are relatively broad because inelastic spin relaxation leads to suppression of the 
^'-matrix poles [54] . 

The magnetic field dependence of the scattering cross sections for collisions of molecules in 
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the absolute ground high-field-seeking state \Ms^ = — 1, M^g = — 1) is dramatically different 
(see Figure 4) . The s-wave elastic scattering cross section displays a manifold of resonance 
peaks which we attribute to the combined action of the interaction potential and the spin- 
spin interaction (1). As a result of an interplay between these interactions, avoided crossings 
occur between the incoming scattering state \Msf^ = —l,Msj^ = ~1) the quasibound 
states of the O2-O2 complex shown in Fig. l{b) leading to the resonant variation of the 
elastic cross section. This suggests that the resonances depicted in Fig. 4 are similar to 
the magnetic Feshbach resonances in coUisions of the alkah metal atoms [12] and that the 
s-wave scattering length in an ultracold gas of molecules can be efficiently manipulated 
with magnetic fields. 

In order to elucidate the possibility for evaporative cooling of molecules, it is useful 
to analyze the dependence of the spin relaxation rates on the rotational constant [20, 55] . 
In Fig. 5, we plot the ratio of the cross sections for elastic scattering and spin relaxation 
as a function of the sphtting between the rotational levels. All other parameters of the 
Hamiltonian were unchanged in this calculation. The probability for spin relaxation in 
atom-molecule and molecule-molecule collisions decreases with increasing Be because the 
matrix clement between the different spin states due to the spin-spin interaction scales as 
Ass /-Be [20, 55]. Figure 5 shows that the elastic to inelastic ratio at a collision energy of 0.1 
K improves dramatically with increasing the rotational constant. The dependence on Bg 
is, however, not monotonic because of the manifold of Feshbach resonances (see Figs. 3-4), 
which are most pronounced at low collision energies. The results shown in Fig. 5 indicate 
that evaporative cooling of molecules with larger rotational constants should generally be 
more efficient. 

Previous theoretical and experimental work has established that the cross sections for 
spin-changing transitions in molecules are sensitive to the splitting between the N — 
and the N — 2 rotational levels. The approximate 1/B'^ scaling law [55] suggests that spin 
relaxation occurs at a slower rate in collisions of atoms and molecules with larger rotational 
constants. This result is important because the sphtting between different rotational levels 
can be altered with electric, microwave, or off-resonant laser fields, thereby enabling control 
over the spin degrees of freedom of cold molecules. We illustrate the idea using the example 
of a nonpolar ^E molecule such as O2 in the presence of far off-resonant, continuous-wave 
laser field. The energy levels are determined by the Hamiltonian (3) with an additional 
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term [57] 

Vl^-j [Aq;cos'x + «±] , (31) 
where Aa = an — q;_|_ is the polarizabihty anisotropy, x is the angle between the molecular 
axis and the laser polarization vector (which we assume to be parallel to B), and Eq is the 
amphtude of the laser field. Figure 6 shows the energy levels of ^^©2 as functions of the 
laser field strength (we transform away the x-independent terms by properly choosing the 
zero of energy) . The effective rotational constant defined as the sphtting between the closest 
field-dressed Zeeman levels in the N = and N = 2 manifolds, increases by a factor of 1.6 as 
the laser field intensity is varied from zero to 5 x 10^^ W/cm^. For polar molecules, similar 
effects can be induced by dc and microwave laser fields [56]. 

As shown in Fig. 6, any substantial modification of molecular energy levels requires 
laser field intensities of order 10^^ W/cm^. The maximum cw laser field intensity currently 
available in the laboratory is about 2.5 x 10^ W/cm^ for a 0.5 mm size sample. However, 
if the trapped cloud is compressed to one-tenth its original size, the maximum attainable 
laser intensity increases by a factor of 100, so the required field strengths of 10^^ W/cm^ are 
well within reach for microscopic clouds (10~^ mm). Samples of such size can be produced 
and stored in magnetic microtraps recently demonstrated for cold polar molecules [60]. 
An alternative solution is to use pulsed lasers, whose intensity is typically much stronger 
than that of cw lasers. However, the time-independent description used in this work can 
only be applied if the duration of the laser pulse is long compared to the collision time. 
This condition is necessary to ensure the validity of the quasistatic approximation for the 
molecule-field interaction [57]. 

Figures 7 and 8 show the cross sections for spin relaxation and elastic scattering in O2- 
O2 collisions as functions of the rotational constant at different collision energies. At the 
highest energy of 0.1 K, the cross section for spin relaxation decreases monotonically and 
that for elastic scattering increases with increasing B^. For a molecule with Bf./Bo^ > 
10, the spin relaxation would be suppressed by two orders of magnitude. As the collision 
energy decreases, the variation of the cross sections with Be becomes more complicated. 
At the lowest collision energy of 10~^ K, both the elastic and spin-changing cross sections 
decrease with increasing Be, and the suppression of spin relaxation in very efficient for 
large rotational constants. In contrast, for the intermediate kinetic energy of 1 mK, the 
probability for spin relaxation remains large even for Bg/ Bq^ > 10, and the cross sections 
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show oscillations due to Feshbach resonances (see Fig. 4). The results shown in Figs. 7 and 8 
indicate that off-resonant laser fields may be used to suppress spin-changing transitions and 
induce Feshbach resonances in ultracold collisions of molecules. We note that it order to 
establish that the laser-induced magnetic Feshbach resonances do occur at practicable laser 
field intensities, it will be necessary to perform accurate quantum scattering calculations in 
superimposed magnetic and off-resonant laser fields. Such calculations can be carried out 
using the formalism presented in this work. 

IV. SUMMARY 

We have presented a quantum mechanical theory and converged quantum mechanical 
calculations for collisions of identical molecules in a magnetic field. The cross sections 
for spin relaxation at magnetic fields above 1 mT are comparable to the elastic cross sections, 
and they are not very sensitive to the applied field. At low magnetic fields, the spin-fiipping 
transitions are suppressed by d-wave centrifugal barriers in the outgoing collision channel, 
suggesting that molecules can be evaporatively cooled in shallow magnetic traps. Our 
calculations shown in Fig. 4 suggest that the cross sections for collisions of molecules in 
the absolute ground state display a dense spectrum of Feshbach resonances, which may be 
used to tune molecule - molecule scattering lengths in a very wide range of magnetic fields. 
These resonances can be used to create tetra-atomic molecules via magnetoassociation [19] , 
manipulate molecule- molecule interactions in optical lattices [7], control chemical reactions 
of polyatomic molecules [4] , or facilitate evaporative cooling of molecules in an optical dipole 
trap [61]. Our results show that the number of Feshbach resonances for molecule-molecule 
scattering in the absolute ground state (~ 100 T^^) is larger by a factor of 10 than the 
density of the resonances for the low-field-seeking states. 
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APPENDIX A: MATRIX ELEMENTS OF THE INTERACTION POTENTIAL 



In principle, the radial coefficients Vxj^Xj^x{R) can be obtained by inverting Eq. (15) 
using the orthonormality properties of the spherical harmonics. This procedure requires the 
evaluation of six-dimensional integrals over the angles {R,rA,rB), which poses a difficult 
practical problem. A more convenient method of evaluating the expansion coefficients in 
Eq. (15) involves a transformation to the body-fixed (BP) frame. The 2;-axis of the BF 
frame coincides with the vector R [40, 41]. The tensor product in Eq. (16) is invariant 
under rotations of the coordinate system [45, 46], so we can substitute R — into Eq. (15) 
and use Eq. (16) to obtain 

1/2 



y,(i2,fA,fB) = (47r)^/^ Y: (^) W(i?) 



X 



^ , Aa Ab A I Y,^^{f'A)Y,,,.m{r'B), (Al) 



m —m 



where the primes indicate the BF angles. In deriving Eq. (Al) we have used the fact that 
Yxruxi^) = [(2A + l)/47r]V2. We define the BF angular basis functions [31, 40, 41] 

Ax.x.xif'A, r'^) = (2A + 1)'/' i ] Y,Uf'A)yx,,-m{r'B), (A2) 

„ ym -m Q J 

to expand the interaction potential as follows 

Ki(i2,rA,rB)=47r ^ AxMir'A,r'B)VxMiR)- (A3) 

Xa,Xb,X 

We note that the radial coefficients V^^Aba(-R) in Eqs. (A3) and (15) are exactly the same. 
This is not the case for the basis functions expressed in different coordinate systems. 

The BF basis functions span a coordinate subspace defined by 4 angles, r'j^ — {9a, <Pa) and 
^B — (^B, </?b)- However, the interaction potential depends on the dihedral angle (f — (fiA — ^PB 
and not on the azimuthal angles separately. By separating the sum in Eq. (A2) into three 
contributions from m > 0, m = 0, and m < 0, respectively, and changing the sign of the 
summation indexes with m < 0, we obtain 



^AAABA(r;,r^)^ ^ J2\ ^ _ \QxMOA)exUOB)[2-dm,o]{-rcos{m^), 

m=0 



m —m 



(A4) 
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where ©AB,-m(^B) are the normahzed associated Legendre polynomials. Here, we have used 
the relation (_)Aa+Ab+a _ ^Yiai, holds for two identical homonuclear molecules. 

In practice, it is more convenient to use the basis functions that depend on the angles 
^A; ^B; and Lp. Even though the right-hand side of Eq. (A4) depends on the three angles, 
the functions on the left-hand side are defined in a four-angle space. Therefore, we need to 
renormalize Eq. (A4) to exclude the integration over one extra angle. We define a reduced 
orthonormal basis set of functions, which span the three-dimensional subspace of angles 

1 /'2A + 1^ 



^E l ^ \®X^m{eK)Qx^m{eB)[2-5m,o]{-rCOs{m^) (A5) 

m=o y rn -m J 

We emphasize that unlike Eq. (A4), this expression defines an orthonormal basis in the 
subspace of three independent angular variables. This suggests the following expansion 

V4R,9a,9b,v)= Y1 VMiR)^XAXnxiOA,9B,^) (A6) 

Aa,Ab,A 

We note that the coefficients VaaAba(-R) are not the same as in Eq. (15). Multiplying Eq. 
(A6) by Ay^x'j^x'i&A, ^b, f) and integrating over all angles, we find 

Kx^xiR)^ deAsm0A d0Bsm0BAx^x,x{dA,OB,^)V{R,0A,eB,<f) (A7) 

Jo Jo Jo 

Comparing the expansion (A6) with Eq. (A3) and using Eq. (A5), we obtain the relation 

between the original coefficients defined by Eqs. (15), (Al), (A3) and those given by Eq. 

(A7) 

VxMiR) = -^KxM- (A8) 



To summarize, the expansion coefficients VaaAba(-R) can be obtained by (i) evaluating the 
reduced basis functions (A5); (ii) integrating the interaction potential with these functions 
following Eq. (A7), and (iii) using Eq. (AS) to obtain the original expansion coefficients 
(15). 
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APPENDIX B: DIFFERENTIAL SCATTERING CROSS SECTIONS FOR IDEN- 
TICAL MOLECULES 



The differential cross section (DCS) for collisions of indistinguishable molecules can be ex- 
pressed in terms of the symmetrized scattering amplitude q^^^^^y ^ {Ru R) given by Eq. (24) 

The integral cross section can be obtained by integrating the DCS over the coordinates of 
the final fiux and averaging over all possible directions of the initial fiux. 



1 



7a7b^7a7b ~ 7Z I I dRidR ^tl" , {Ri, R) (Bl) 



7. 

4^yy ---i--- ^^^^^ 

Since the scattered molecules in the same internal state are not distinguishable after the 
collision, the integration over R in Eq. (Bl) has to be restricted over half-space [25, 49] for 
the final states satisfying 7^ = 7b- Thus, Eq. (Bl) can be written as 



1 



27r pw /•2n /-^max da 



<^7A7B-7;7^, = 1^ / d<Pi I sm9id9i I d(p I sinOde ""'"f (^i, 0i; 0), (B2) 
^ ^ 4:n Jo Jq Jq Jq dRidR 

where Ri = (6'i,0i), R = {0,(f)), and we have defined ^^^ax to be 7r/2 if 7^ = 73 and tt 
otherwise. The integration in Eq. (B2) can be performed trivially for ^max = tt to yield 
[40, 41] 



^7A7B->7i7^5 - ^"T^ -^^^ \^lKl^tme,lWm'^ 1^ 7^ 7b- (^3) 

7A7B i^rai I'm'^ 

All that remains is to consider the special case 7a = 7b- Expanding the square modulus in 
Eq. (B2) yields 



_ 7r(l + V^b)2 

<^7a7b^7a7b U2 

7A7B 



^ ^ ^ ^ j-^i+^'i+^a-^^ 
■ti,me^ £2,m(^ t\,m'^^ ^2 '"^^2 



27r 



^^'l^B4m.,;7Mm^7'-^B^^-^2;7i7^,^w^J^ ^'^i I sin^id^ir,,^,^(^i,0i)5^,;„^,(^i,0i) 

The first integral is the usual normalization integral of two spherical harmonics. The second 
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integral can be separated in two parts [45] 

f27r /■7r/2 /■27r 



X sin ^ci^e,.„.^ (^)©^W^^ (^)- (B5) 
The one-dimensional integral over (p is equal to S^' m' , so that the intergal (B5) becomes 

m', r' sinedeee^m', (^)e,^„^ (^). (B6) 

^ ^ Jo ^ ^ 



Since for identical molecules and £'2 are of the same parity (either even or odd) [62], and 
the index m^^ is fixed, the associated Legendre polynomials in Eq. (B6) are both even or 
odd functions of cos^ [45, 46]. Therefore, their product is always an even function of cos^. 
Since in this case the substitution 9 ^ n — 9 leaves the integrand in Eq. (B6) unchanged, 
we can halve the 9 integration range and obtain 

d0 sin 9d9Ye,^,^ {9, ct>)Yl^,^^ {9, <t>) = \5e,eM,rK, ■ (^7) 

Substituting this result into Eq. (B4) and performing the integration leads to the cancella- 
tion of the prefactor (1 -|- (^7^7^) by the | from Eq. (B7). This completes the derivation of 
Eq. (29). 
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Figure 1: (a) Energy levels of the ^''02(^S^) molecule as functions of the magnetic field, (b) A 
schematic illustration of spin relaxation in collisions of ^T, molecules. 
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10"** 10"'' lO'" 10"^ 

Collision energy (K) 

Figure 2: The cross sections for elastic scattering (a) and spin relaxation (6) in 02('^S~) - 02(^5]~) 
collisions as functions of the collision energy at different magnetic fields: 10~^ T (circles), 10~^ 
T (squares), and 0.1 T (diamonds). The ratio of the cross sections for elastic scattering and spin 
relaxation is shown in panel (c). 
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Figure 3: The cross sections for elastic scattering (full line) and spin relaxation (dashed line) in 
collisions of 02(^S~) molecules in the low field-seeking state \Msj^ = l,Msg = 1) as functions 
of the magnetic field. The spin relaxation cross section is summed over all final spin states and 
divided by 10 to fit the scale of the figure. The collision energy is 10^^ K. 
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Figure 4: The s-wave elastic scattering cross section for collisions of ©2(^5] ) molecules in the 



lowest high- field-seeking state \Msj^ 
collision energy is 10~^ K. 



-1) as a function of the magnetic field. The 
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Figure 5: The ratio of the cross sections for elastic scattering and spin relaxation in collisions of 
02('^5]^) molecules in the low-field-seeking state jM^^ = 1, Ms^ = 1) as a function of the rotational 
constant at different collision energies. The magnetic field is 1 T. 
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Figure 6: Energy levels of ^^02(^S~) in superimposed magnetic and off-resonant laser fields. The 
levels are shown as functions of the laser field intensity at a fixed magnetic field of 1 T. The dashed 
line shows the highest low-field-seeking state \Msj^) = 1. The laser field intensity is defined as 
Iq = Eq (see text for details). 
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Figure 7: The cross sections for spin relaxation in collisions of 02(^S~) molecules in the low-field- 
seeking state \Ms^ = l,Ms^ = 1) as functions of the rotational constant at different collision 
energies. The magnetic field is 1 T. 
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Figure 8: The cross sections for elastic scattering in collisions of 02(^S^) molecules in the low- 
field-seeking state \Msf^ = l^Ms-^ = 1) as functions of the rotational constant at different collision 
energies. The magnetic field is 1 T. 
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